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on : Recent observations have revealed that small-scale structure exists at high redshift. We 

in . 

study the possibility of the formation of such structure during baryogenesis and big bang 
nucleosynthesis. It is known that under certain conditions, high density baryonic bubbles are 

in ■ 

created in the Affleck-Dine model of baryogenesis, and these bubbles may occupy a relatively 
small fraction of space, while the dominant part of the cosmological volume is characterized 



o 



by the normal observed baryon-to-photon ratio, r\ = 6 ■ 10 -10 . The value of r] in the bubbles 
could be much larger than the usually accepted value (indeed, even close to unity) and still 
be consistent with the existing data on light element abundances and the observed angular 
spectrum of cosmic microwave background radiation (CMBR). We find upper bounds on 
r) by comparing the abundances of heavy elements produced in big bang nucleosynthesis 
(BBN) and those of metal poor stars. We conclude that 77 should be smaller than 10~ 5 in 
some metal poor star regions. 



X 



§1. Introduction 

One of the biggest puzzles in astrophysics is the nature of structure formation 
from early stages. We know that there exist galaxies, quasars (QSO) and stars in 
our universe. These structures must have formed at some stage in the evolution 
of the universe. But when did the universe become light? Recent observations 
suggest that such structure formations began earlier than conventionally thought. 
For example, Wilkinson Microwave Anisotropy Probe (WMAP) data suggest that 
reionization began when z ~20. 1 ) Also according to Rets. 2 ), 3 ' star formation activity 
started when z ~ 10, or at some slightly large value. In addition, it is known that 
the quasar metallicity did not significantly change from the time of high redshift to 
the present time. 4 ) Some quasars already reached solar or higher metallicity when 
z was no smaller than 6. Recently a galaxy at £=10.0 was observed. 6 ) For IGM, 
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the bulk of metal ejection occurred when z > 3. 7 ^ With regard to the abundances of 
metal poor stars , see Ref. 8 ) 

Do these QSO, intergalactic medium (IGM) and stars from high redshift con- 
tradict the theory of structure formation? 

We still do not have a "standard theory" of stars and QSO formation. But 
in some models (for example, 9 ) 10 )), these observations do not seem to yield strong 
contradictions in the cases of QSO and reionization, while the galaxy observed at 
z=10 corresponds to the collapse of > 2a fluctuations. 6 ) 

We can adjust the theory by changing the initial mass function and other pa- 
rameters in order to account for the observations. However, because we understood 
the mechanism of star formation only poorly and also because observations suggest 
some possibility of the existence of structures at higher redshift, it is valuable to 
consider an alternative scenario for objects in the early universe that have been or 
will be observed. 

We propose a scenario in which the seeds of these structures are produced during 
baryogenesis. One of the present authors found 11 ) that under certain conditions, high 
density baryonic bubbles with small spacial size are produced, while most of the 
universe is characterized by the small baryon-to-photon ratio (77) observed through 
BBN 12 ) and CMBR, 13 ) rj = 6 • 1CP 10 . The model represents a modified version of 
the Affleck-Dine 14 ) baryogenesis scenario and is based on the hypothesis that the 
Affleck-Dine field (j) is coupled to the inflaton field The interaction Lagrangian is 
assumed to have the general renormalizable form 

£ int = \\<P\ 2 <P 2 + g\<P\ 2 4> 

= A(<2>-<2> 1 ) 2 |</>| 2 -A<2> 2 H 2 , 

where g and A are the coupling constants and <Pi = —g/2X. It is known that the 
effective mass of the field 4> may contain the contributions 

( m f 7/ ) 2 = ml + £R + f3T 2 + A(# - <2>!) 2 , (1-2) 

where £ and (3 are constants, R is the curvature scalar and T is the temperature of the 
primeval plasma, rrig is the vacuum mass of (j), barring the contribution from (Ai^i) 2 . 
For minimal coupling of (ft to the gravity field, the coupling to the curvature vanishes, 
by definition. However, radiative corrections may induce relatively small coupling, 
with £ < 10 -2 . Temperature corrections to the mass appear at higher orders in the 
perturbation theory, and thus usually <C 1. In what follows, we ignore these two 
contributions, as they are not essential for the dynamics of the formation of high 
B bubbles (but may be important for some details of their evolution). We assume, 
though it is not necessary, that m\ ~ H 2 , where Hj is the Hubble parameter during 
inflation. The coupling constant A is bounded from above by the condition that 
the inflation self-interaction be sufficiently weak. The interaction (1.1) induces the 
inflaton self-coupling ~ A 2 ^ 4 , and the condition A ~ 10 -5 -10 -6 allows it to remain on 
a safe level to avoid density perturbations that are too large. An essential condition 
for the realization of the high-B-bubble scenario is a negative effective mass squared 
of the field 4> when the inflaton field, 0, is close to <P\. To this end, we stipulate the 
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condition 

mg+£R<0, (1-3) 

and because we assumed that the term £R is negligible, this condition is satisfied 
for negative TOq. In fact, a weaker condition is possible, namely that < Hj(<P\), 
where Hi{^>i) is the Hubble parameter for ^ ~ (Pi, and it is assumed that inflation 
still proceeds at this stage. This condition is necessary, in particular, to ensure a 
sufficiently large size of the bubbles. 

As argued in Ref. 11 ) , the scenario of the formation of bubbles with high baryon 
asymmetry would be relevant, for example, for Hj ~ (10 -5 -10 -6 )mp;, A ~ 1CP 10 - 
10~ 12 , and <Pi ~ mpi, but other wider ranges of parameter values also seem to be 
allowed. 

The inflaton field <P is assumed to evolve from some high value 3> > mpi down to 
zero, where inflation ends. In the course of this process, the effective mass squared 
of (p changes from a positive value to a negative value (or to \m e ff \ < Hj), and then 
back to a positive value. 

An important assumption is that the potential of the field possesses two min- 
ima at small values of m e ff. This is true, e.g., for the Coleman- Weinberg 15 ) potential 

U{4>)=ml }} \ ( p\ 2 + \^\ i \n ] ^. (1-4) 

We assume that chaotic inflation is valid, though not a necessary condition. Initially 
when the inflaton field has a large value, i.e. <P > mpi, and the square of the effective 
mass is positive, the potential (|1-4I) has only one minimum, at <p = 0. As <P decreases, 
the effective mass decreases, and a new minimum at <j> ~ a appears. When 
becomes negative or zero, the minimum at = becomes a local maximum, and the 
field (j) tends to increase from to a larger value in the direction of the minimum at 
If m e ff remained negative for a sufficiently long time, the field cf> would tend 
to this other minimum in the entire space. However, because m^ff remains negative 
for only a finite time, depending upon the duration of the "negative" period and the 
magnitude of the fluctuations of (f) near (ft = (the latter being typically of order 
5(ft ~ H), only some bubbles with large values (j) ~ °~ occupying a small fraction of 
space are formed, while the rest of space is characterized by <p = 0. 

As <P decreases further, the effective mass becomes positive again, and the second 
minimum at 4> ~ a becomes higher than the first one at <p = 0, and ultimately it 
disappears. Subsequently the field <p stuck in this minimum decreases to zero. If this 
field carries baryonic charge, as in the Affleck-Dine scenario of baryogenesis, then 
a large baryon asymmetry is generated inside such bubbles, while in the region of 
space occupied by eft near zero, the asymmetry is much smaller. This asymmetry 
may be generated by the decay of the same <f> field as in the bubbles but with a 
smaller amplitude or by some other mechanism of baryogenesis, which is usually less 
efficient and creates an asymmetry much smaller than that in the Affleck-Dine case. 
It is worth noting that the scenario of baryognesis inside the bubbles suggested in 
Ref. 11 ) and used here is very similar to the original Affleck-Dine scenario with the 
only difference being that in the original version, the field (j) is displaced from the 
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minimum at <f> = along the flat directions of the potential by quantum fluctuations 
during inflation, while here the field <f> is displaced from zero to a large value by a 
negative m^ff- Here it is worth repeating that because rnlff rema i ns negative during 
only a short time interval, only a small fraction of is able to reach the minimum 
near the large <p = a. 

As a result of the process described above, the bulk of the universe would come 
to possess the normal small baryon asymmetry observed through BBN and CMBR, 
7] = 6 ■ 10 -10 , while in some small regions, the asymmetry could be much larger, 
perhaps even reaching r\ ~ 1. Depending upon the details of the 4> potential and 
the mechanism of CP-violation (though in some versions of the scenario, explicit 
CP-violation is unnecessary), the value of the asymmetry may vary form bubble to 
bubble or remain the same, in which case only the bubble size may vary. The baryon 
asymmetry inside the bubbles cannot be reliably predicted because of the existence 
of many unknowns in the theory, but typically, it should be the same as in the usual 
Affleck-Dine model. 

After baryogenesis, the initial energy density difference between the interior of 
the bubbles and external space is small (isocurvature perturbations). When the 
QCD phase transition takes place and quarks form nonrelativistic baryons, large 
density inhomogeneities develop, because the equation of state inside the bubbles 
begin to deviate from the relativistic one, p = p/3, which is valid in the rest of the 
universe. The subsequent destiny of these high density baryonic bubbles depends 
on the size of the bubbles and the value of the baryonic charge asymmetry. Some 
of them may form unusual stars or anti-stars with a high initial fraction of heavy 
nuclei, and others may form primordial black holes. 

According to the calculations of Ref. 11 ) , the mass distribution in these region 

is 

dN ( ~, 2 M\ 

dM (XeXP (" Cln Mi J' (1 " 5) 
where C and M\ are constants which should be determined by Hj, the time evolution 
of m e ff when <P ~ <Pi, the time width <P = <&\ and <P = 0, etc., which in turn are 
determined by the unknown details and parameters of the potentials of the inflaton 
and <fi fields. "Natural" values of these constants are C < 1 and M\ = (10 -3 - 
1O 6 )M . The distribution of bubble sizes is similar. We assume that the bubble 
sizes are smaller than the galaxy size. 

Because this model has a great degree of freedom and the parameters can be 
treated as free, we cannot a priori restrict such important quantities as, e.g., the 
bubble size and the baryon density inside the bubbles. Instead, assuming that this 
model is applicable to the early universe and studying its implications for the for- 
mation of primordial objects, dark matter, elemental abundances, etc., one could 
either find observational confirmation of the discussed mechanism in the generation 
of large baryonic inhomogeneities at small scales or obtain bounds on the magnitude 
of the effects and thus on the parameter values used in the model. 

In order to carry out either of the above stated tasks, we need observational 
data concerning BBN and CMBR. Observations of the abundances of primordially 
produced light elements 12 ) allow the "measurement" of the baryon-to-photon ratio 
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V = ns/n 7 during BBN. On the other hand, the spectrum of angular fluctuations 
of CMBR 13 ) also yields a value of 77 that is in good agreement with BBN result. 16 ) 
Based on the observational BBN and CMBR data, many interesting restrictions on 
models of baryogenesis, the magnitude of the reheating temperature, unobserved 
particle species, their masses, lifetimes, and coupling strengths, possible types of 
phase transitions in the early universe, etc., have been obtained. Most of these 
studies deal with standard (or homogeneous) big bang nucleosynthesis (SBBN), 17 ) or 
inhomogeneous big bang nucleosynthesis (IBBN), 18 -* but with a very small magnitude 
of the fluctuations of the baryonic charge density. However our case is different. In 
our model, large inhomogeneities in spatially small regions are essential, and from 
the observation of small scale structure, we determine restrictions on the nature of 
elementary physics. Though CMBR supports the existence of scale invariance in the 
primordial power spectrum, recent observations suggest it is reasonable that there 
was small scale and small fraction of structure in the early universe. Motivated 
by these findings, we consider here the hypothesis that the presently existing metals 
could have been produced in the very early universe during BBN. We have computed 
the abundances of different (not only light) elements produced in the bubbles with 
high baryon density and compared the abundances of heavy elements produced in 
these bubbles during BBN with those in QSO, IGM and metal poor stars observed 
at high redshifts. 

§2. Heavy element production at BBN with a large baryon number 

We assume that the characteristic bubble size is much larger than the baryon 
diffusion length and hence that baryon diffusion is not important. In this case, the 
problem is greatly simplified, because we need only consider homogeneous nucleosyn- 
thesis. The reaction network used in this work was applied previously to supernova 
nucleosynthesis calculations 19 ) 20 ) , but we add 16 light nuclei. The network includes 
the 258 isotopes listed in Table 1, Z = - 32 and A = 1 - 74. 

The initial and final temperatures are 10 11 K and 10 7 K, which correspond to the 
time interval from 10 -2 to ~ 10 6 sec. 

In our calculation of the time evolution of the baryon density and temperature, 
we use the Friedmann equation 

H 2 = ^Gptoua, (2-1) 
where ptotai = P-y + (p e - + Pe+) + Pu + Pb-, and the energy conservation law 

where the last term takes into account the change of energy introduced by nucle- 
osynthesis. We do not consider the possibility of neutrino degeneracy. 

Calculations of nuclei production during BBN with high rj have already been 
carried out by one of the present authors. 21 ) There are basically two distinctions 
between the present work and the that work. The first is that the network used in the 
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Mn54 Fe55 Co55 
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Fe60 C06O Ni60 Cu61 

Fe61 C06I Ni61 Cu62 

Co62 Ni62 Cu63 

Co63 Ni63 Cu64 

Co64 Ni64 Cu65 

Ni65 Cu66 
Cu67 
Cu68 



Table I. Isotopes that are included in our nuclear reaction network 
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previous work includes 72 isotopes, while ours includes 258 nuclei. Therefore, we can 
predict which heavy elements should be observed, unlike in the case of the previous 
work. The second is that the main interest of the previous work was to study the 
effects on BBN of a large cosmological lepton asymmetry with a magnitude greater 
than that of the baryonic asymmetry, while we assume that the lepton asymmetry 
is negligibly small. 

§3. Results and Discussion 

In our model, the only free parameter is rj. We carried out the calculation for 
values of t] ranging from 1CT 12 to 3 x 1CT 4 and. Table 2 displays the mass fractions 
of the main product elements. 
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Table II. Main product nuclei and their mass fractions for r\ = 1CT 4 , l(r 6 and 10" 10 . This result 
shows that the amount of heavy elements produced increases as ij increases. We also see that 
many of the product nuclei are proton rich. 

As 7] becomes larger, heavy elements begin to be produced more efficiently. A 
very interesting feature of our result is that the nuclear reactions proceed along the 
proton rich side. The usual nucleosynthesis of heavy elements in supernovae proceeds 
along the neutron rich side (r-process). However our calculations show that BBN 
produces proton rich nuclei. Though we cannot conclude that BBN is a p-process, 
we can say that most of the product elements are proton rich. To confirm this, we 
calculated the value of Ye= n p /(n p + n n ). The results are presented in Table 3. 

It is apparent from Table 3 that the nuclear reactions proceed in a proton rich 
environment. The value of Ye is determined by the amount of He 4 produced, and 
hence the (3 decay effect is almost negligible. The r\ dependence of each element 
abundance is displayed in Fig. 1. 

The abundances of H and He remain almost constant over a wide range of values 
of r). C and O increase as 7] increases, but eventually they reach maximum values 
and then decrease. Beyond that point, heavier elements, such as Ca and Fe, begin 
to dominate. Two typical examples of the time evolution of the abundance of each 
element are plotted in Figs. 2 and 3. 

In the early stages, there is an abundance of light elements. Then as reactions 
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Table III. Ye for r/ — 10 10 - 10 5 . It is clear that the nuclear reactions proceed in a very proton 
rich environment. This is the reason for the abundance of proton rich nuclei. 
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Fig. 1. 77 dependence of the mass fraction of each element. These of H and He are almost constant. 
Contrastingly, heavy elements are produced in greater abundance as rj increases, but eventually 
they begin to be consumed to produce even heavier elements. 
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Fig. 2. Time evolution with 77 = 10 
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time(9ec) 

Fig. 3. Time evolution of nuclei abundances with r\ — 10~ 8 

proceed the light elements are consumed, and heavy elements begin to be produced. 
Production of heavy elements at rj = 10 -5 is much larger than at n = 1CP 8 . In 
particular, nuclei with mass numbers greater than 73 begin to dominate at rj = 10~ 5 . 

In Fig. 4 the ratio of the BBN number fraction with respect to the present 
solar system number fraction is presented. We find two peaks in the graph, one 
corresponding to calcium at Z=20, one corresponding to Germanium at Z=32. We 
believe that the peak value at Z=32 should not regarded as the Ge fraction but, 
rather, as the sum of the abundances of all elements with atomic number greater 
than 32, because our network includes only those nuclei with Z<32, and thus the 
reaction cannot proceed beyond Germanium, and the nuclei accumulate there. By 
contrast, the peak at Z=20 should be considered as representing a peak of calculation. 
To understand why Ca is produced to much a great degree, it is necessary to know 
which of the Ca isotopes is produced in greatest abundance. These data are presented 
in Table 4. The table shows that the dominant output is for Ca 40, which has 20 



isotopes abundances 



Ca40 


5.33 x 10" 


-8 


Ca41 


1.8 x 10 _ 


9 


Ca42 


1.93 x 10- 


-9 


Ca43 


1.21 x 10" 


-12 


Ca44 


1.89 x 10" 


-17 



Table IV. The abundances of Ca isotopes. This shows that Ca 40, which is the double magic 
number nucleus, is dominantly produced. 

protons and 20 neutrons. It is to be noted that 20 is a magic number, and Ca 40 is a 
double magic number nucleus which is known to be very stable. Surprisingly, no peak 
can be found at Z=26 (Fe). On the contrary, there is a decrease there. There are 
two reasons for this. The first is that because of the p-process, the production of Fe 
is inhibited. The second is that the solar abundance of Fe is large. At rj = 3 x 10 , 
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[Fe/H]= —6.70 is still too small, but [Ca/H]= —2.17 is large enough to be inconsistent 
with the observations of extremely metal poor stars. For example, [Ca/H]= —5.37 
for He0107-5240. 22 ) In order for the Ca abundance produced at BBN to be below 
this value, rj of a high baryon density bubble, if surrounding this metal poor star, 
should be smaller than 10 -5 . However, there is such a great degree of freedom in our 
model that we cannot strongly restrict it with data from only a single observation. 
We can only say that if such a model is realized, rj of the bubble around this metal 
poor star must be under 10 -5 . On the other hand, it may be the case that the 
observed early metal poor stars are well outside of the baryon rich regions. For 
confirmation of the model, we need to carry out calculations with larger values of 
rj, because for rj < 3 x 10~ 4 , the abundances of heavy elements are still small in 
comparison with those observed in IGM, around QSO, and in most of the metal 
poor stars. Observations of objects in the early universe with anomalously high 
metal abundances are also desirable. To impose some restrictions on the mass of 
the Affleck-Dine field, its coupling strength, etc., we need additional observational 
data. Extending our calculations to larger values of rj in order to obtain predictions 
for the abundances of heavier elements, it may be possible to distinguish usual stars 
from our bubble-made stars, which is necessary to improve the restriction. With 
such progress, in the future, we may be able to reach a more definite conclusion. 

There is another interesting possibility. To this point, we have considered 
the case in which baryons in the bubbles do not diffuse. However, if the bub- 
bles are smaller than the quark diffusion distance, which can be evaluated 11 ) in 
comoving coordinates as ~ (tlfree) 1 ^ 2 : where lf ree is the quark mean free path, 
I free ~ {°-N)~ l ~ T -1 , they do diffuse. This could affect the angular spectrum of 
CMBR. Most IBBN studies carried out to this time treat only small values of i], 
specifically, r] ~ 10~ 10 . However, with our approach, we are able to investigate a 
novel relation between IBBN and CMBR. Our study also suggests some possibili- 
ties for the origin of p-nuclei in the solar system. To investigate this, we need to 
calculate the BBN abundances with a reaction network that includes heavier nuclei. 
There is also the possibility that the baryon rich bubbles, though not forming pri- 
mordial black holes, might end up as stellar-type or planetary-type objects. In this 
case (which is under investigation), the observational consequences could be quite 
different. 
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